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We consider the evolution of narrow-band wave trains of 
finite amplitude in a nonlinear dispersive system which is 
described by the Klein-Gordon equation with arbitrary polynomial 
nonlinearity. We use a new perturbative technique which allows the 
original wave equation to be reduced to a model equation for the 
wave train envelope (high-order nonlinear Schrodinger equation). 
The time derivative is expanded into an asymptotic series in 
two independent parameters which characterize the smallness of 
amplitudes (e) and the slowness of their spatial variations (^t). 
In contrast to other perturbative methods in which these two 
parameters are taken equal (e.g., the multiple scale method), the 
two-parameter method produces no secular terms. The results 
of this study can be applied to investigating the propagation of 
ultrashort (femtosecond) pulses in optical fibers, to studying the 
wave events on a fluid surface, and to describing the Langmuir 
waves in hot plasmas. 



1. Introduction 

One of the remarkable properties of nonlinear dispersive 
systems is the possibility of the existence of steady 
progressive waves with finite amplitude due to the 
balance between dispersion and nonlinearity effects. 
In the absence of dissipation, the dispersion relation 
(dependence of wave frequency on wavelength) can be 
written as 



u=u{k, \A\ 2 ), 



(1) 



where k = 2ir/\ is the wave number, uo is the wave 
frequency, and A is the wavelength. In systems with no 
dispersion, the dependence on k is linear, i.e., u> = ck. 
The constant c is called the wave phase speed. In the case 
of nonzero dispersion, the phase speed depends on wave 
number (d 2 uj/dk 2 ^ 0), so that wave trains spread out 
in space. In nonlinear systems, this dispersive spreading 
can be compensated by nonlinear effects, which manifest 
themselves in the dependence of wave frequency on wave 
amplitude A in dispersion relation |T]). 

For narrow-band wave trains, when Afc <C k, the 
dispersion relation can be expanded into a Taylor series 



about the carrier wave number fco and frequency ujq: 

1 fd 2 > 

,k + 

k—kn 



Ak 2 - 

k—kn 



( dto 



\A\< 
=0 



+ 



Going over to the operator equation for the amplitude 
A by the substitution Acu id/dr, Ak —id/dS, and 
omitting the high-order terms, we get 

i(A T + aiAt) - a 2 A^ + a , , A\A\ 2 = 0, (3) 

where a x = (du/dk) k=ko , a 2 = \{d 2 u) / dk 2 ) k=ko , 
ao,o,o = (<9 W /5|^| 2 )|A|=0) an d the variables r and £ 
stand for some conventional time and coordinate with 
respect to which the amplitude of wave train envelope 
exhibits slow variations. 

Equation J3]) is called the nonlinear Schrodinger 
equation (NLSE). It is met, in particular, in various 
problems of nonlinear optics, plasma physics, and 
hydrodynamics, and it admits solutions in the form of 
solitons [1, Sec. 8]. NLSE takes into account the second- 
order dispersion effects (term A^) and the phase self- 
modulation (term ^4|v4| 2 ). These terms are sufficient to 
describe the propagation of picosecond pulses in optical 
fibers [2], Since solitons are formed by the balance of 
dispersion and nonlinearity, a lesser pulse width can be 
obtained when the carrier wavelength is chosen such that 
the coefficient a,2 of the dispersion term vanishes. So, 
high-order dispersive and nonlinear effects come to the 
forefront in femtosecond-pulse problems. These effects 
are described by generalized NLSEs (high-order NLSEs) 
[2,3]. In the case of hydrodynamic wave propagation 
along a fluid surface, high-order NLSEs were derived, in 
particular, in [4-7]. Another distinctive feature of wave 
trains described by NLSE is their instability with respect 
to long-wave modulations (modulational instability) at 
0200,0,0 < [8, p. 640]. To determine the modulational 
instability conditions at ao,u,o = 0, high-order terms are 
to be taken into account in Eq. ([3]). 
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There are several methods which allow the original 
equations of one or other physical process to be 
reduced to model evolution equations of NLSE type. 
They include the multiple scale method [1, 6, 9], 
Hamiltonian formalism (Zakharov method) [10-12], 
variational method [13], and reductive perturbation 
methods [14, 15]. In the multiple scale method, an 
unknown function u(x, t) of coordinate and time is 
looked for in the form of asymptotic expansion in powers 
of a small nonlinearity parameter e: 



u(x, t) 



,(n) 



(x, t). 



(4) 



The wave motion is classified into slow one and fast one 
by introducing different time scales and different spatial 
scales: 



± n — M ? ^ n 



[I X. 



The derivatives with respect to time and coordinate are 
expanded into the following series: 



d_ 

dt 



d 



n— 



d_ 

dx 



E^ 1 



d 



(5) 



the times T n and coordinates X n being assumed to 
be independent variables. A principal drawback of this 
method lies in the fact that the parameters e and \x with 
different physical meanings (the former characterizing 
the smallness of nonlinearity, and the latter describing 
the slowness of temporal and spatial variations) are 
tentatively taken equal: e = //. This admission produces 
the so-called secular terms in the equations for u^ n ' (x, t). 
Such terms, which infinitely grow with time, are 
eliminated in each new order of e by an appropriate 
choice of free parameters emerging in the solutions of the 
linear inhomogeneous wave equations derived from the 
original nonlinear equations for the function u(x, t). The 
procedure is very awkward, and it is difficult to formulate 
in algorithmic form. In Zakharov's method, the problem 
is reduced to an integral equation in the Fourier space, 
and the corresponding solutions should be transformed 
back to the physical space with the use of the inverse 
Fourier transformation. Again, the procedure is quite 
laborious. The same remarks can be made regarding all 
other methods mentioned above. 

In work [16], V.P. Lukomsky proposed an idea of 
constructing a perturbation procedure free of secular 
terms. It was used to derive a generalized NLSE for the 
modulations of gravity waves on deep water. The method 
allows the original system of nonlinear equations to be 



reduced to a model equation for the pulse envelope in the 
form of asymptotic expansion of the time derivative in 
terms of two independent parameters which characterize 
the smallness of amplitudes (e) and the slowness of 
their spatial variations (/j,). In this paper, we present 
a general realization of this two-parameter procedure 
by the example of the reduction of the Klein-Gordon 
equation to a generalized NLSE. Our technique allows 
the coefficients of the generalized NLSE to be calculated 
in arbitrary order of s (high-order nonlinear terms) and 
(high-order dispersive terms) as well as any their 
combination (nonlinear-dispersive terms). 

Consider some wave process described by the (1 + 1) 
Klein-Gordon equation with arbitrary polynomial 
nonlinearity: 



u tt - c 2 u xx + a p u p = 0. 
p=i 



(6) 



Here u is an unknown twice differentiable function 
of the wave process, < t < oo is time, — oo < 
x < oo is coordinate, c and a p are arbitrary real 
constants (a± ^ 0), and P is an arbitrary positive 
integer. Let the initial condition at t = have the 
form u (x , 0) = Q(x) (exp(ifcir) + exp(— ikx)j , Ut(x, 0) = 
P(a;)(exp(ifcx)+exp(— ikx)) , where k is the carrier wave 
number. 

The Klein-Gordon equation arises in the field theory, 
elementary particle physics, crystal dislocation models, 
etc. [1]. When a 2p +i = (-l) p /(2p+l)!, a 2p = 0, P = oo, 
Eq. J6]) is called the sin-Gordon equation, and it is used 
to model the dynamics of dislocations in crystals, self- 
induced transparency in nonlinear optics, spin waves in 
fluid helium, propagation of fluxons in long Josephson 
(superconductive) junctions, and dynamics of domain 
walls in ferromagnetics [8, p. 840]. 

2. Spectral Representation 

We look for a solution to Eq. (J6J in the form of truncated 
Fourier series with variable coefficients: 



N u 



i(x, t) = ^2 Unix, t)e x 



n(u)t—kx) 



(7) 



--N u 



where cu is the wave-train carrier frequency, N u + 1 is 
the number of harmonics taken into consideration, and 
* stands for complex conjugate. The same series can be 
written for all integer powers of the function u: 



vP(x,t)= Y, K)«0M j " 

n=—pN u 



in(iot—kx) 



, P = 2, P, (8) 
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where (it p )_„ = (!t p )*. The coefficients (u p ) n can be 
expressed recurrently in terms of the coefficients u n [17, 
p. 30]: 

min(iV„, n+(p— l)N u ) 
ni=max(-N u , n — (p— 1)JV„) 

The corresponding expansions of the derivatives are 
u tt (x, t) = ^ ((u„)tt + 2incj(u n ) t - 

n=-N„ 



(9) 



BX (a;, t ) = ((ti n ) ra - 2in/c(i 



-n 2 fc 2 u n )e ln( ' , '* _fc3,) . (10) 

Substituting (j7|)- ([T0"l) in © and equating the coefficients 
at the like powers of the exponent exp(i(wi — kx)), we 
obtain a system of nonlinear differential equations for 
the coefficients u n {x, t) (n = 0, N u ): 

{u n )tt - c 2 {u n ) xx + 2in(ui(u n ) t + c 2 k(u n ) x ) + 

p 

+ (n 2 c 2 k 2 - n 2 uj 2 + ai)u„ + «pK)« = 0- (11) 

p=2 

Linearization of these equations at n = 1 gives the 
dispersion relation in the linear approximation: 



2 i 2/2 

oj = oti + c k . 

3. Two-Parameter Expansions for 
Narrow-Band Wave Trains 



(12) 



Generally, the system of equations (fTTj) is by no means 
more simple than original equation It can be 

simplified if solutions are looked for in a class of 
functions with narrow spectrum, |Afc| <§; k (quasi- 
monochromaticity condition). In this case, the problem 
has a formal small parameter fi ~ \Ak\/k, and the 
coefficients u n (x, t) can be regarded as slow functions 
of x and t. Let us introduce a slow coordinate £ = [ix 
and go over to the variables u n — u n (^x, t). 

When there are no resonances between higher 
harmonics, the amplitudes of Fourier coefficients 
decrease with increasing number (quasi-harmonicity 
condition) : 



where u\ = eA. The parameter e can be chosen 
as a second formal parameter, which is independent 
of the dispersion parameter /i in the general case. 
The use of two independent formal parameters is a 
distinctive feature of our approach as compared to other 
perturbative methods (e.g., multiple scale method), 
where these parameters are not distinguished (e = fi). 
When these incomparable parameters are set equal, a 
perturbative procedure produces non-physical secular 
terms. 

In contrast to perturbative methods which use the 
expansions of form |(4]) and ((5j) to reduce Eqs. (fTTjl to 
evolution equations of NLSE type ([3]), we immediately 
start from the most general explicit form of such an 
evolution equation. To this end, the time derivative 
(u\) t = eA t should be expressed in terms of the 
derivatives (ui) nx = e[i n A n £ with respect to coordinate 
(designation A n ^ means the n-th derivative with respect 
to £) and all possible combinations of nonlinear terms 
e 2n+1 A(™ +1 )(A*)". Hence, the derivative A t can be 
written as the following asymptotic expansion in terms 
of parameters e and /j,: 



oo 

A t = iJ2 M 

n =0 



y ] y ] a n -ni, n 1 ~n 2 ., n 2 ^-(n -ni)4 X 
m— ri2— 



(14) 



Expression (fT4"| is the general form of the evolution 
equation for the complex amplitude A of the first 
harmonic. The unknown coefficients a n _ can be 
determined from Eqs. (fTTjl . To this end, the amplitudes 
of all other harmonics (ug, U2, 113, . . .) are expanded in 
terms of the amplitude of the first harmonic A in the 
same manner as it is done in expansion (fl4|) : 



uo 



OO Tl 



no— m= 



no ni no 

V V V h {0) 

/ ; / ; / ; "no— ni.ni- 
ni=0 U2—0 713— 



™2, "2—13, "3 ' 



e n A, n > 1, u ~ e 2 A, e < 1, 



(13) 



xA {no _ ni)i A {ni _ n2H Al n2 _ n3)£ Al 3i + 0(e 4 )) , (15) 

00 no 
U 2 = e 2 (iM)"° ( E 6 £,-m,ni Ano-nOcAiifi + 



no=0 ni— 
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n m n 2 



. f 2 r w h {2) 

/ j / j / j u n -n 1 ,n 1 -n 2 ,n 
ni=0 n 2 =0 ns=0 



2— n 3 , ™3 ' 



xA (n -nx)i A {nx-n 2 )S A (n2-n 3 )i^-n3i + O^ f) > ( 16 ) 
oo no Hi 



"3 

no— ni- 0n2— 

X ^(no-ni)C J 4(ni-n 2 )C^n 2 ? + 
no til n 2 n 3 

+£ 2 y- y- y b (3) 

Til— 712—0 ri3— Ti4=0 



2 , 712 — ™3 ri3 — 714, ri4 ' 



X ^(no-7i.i)^^(7ii-7i2)C^tn2-n 3 )C^(7i3-7i4)C^n4C + ^( £4 )) ' 



(17) 



(n) 

The unknown coefficients bh, are found along with 
the coefficients a n _ from the system of equations l|lip 
by substituting expressions lfl4|) - (fl7|l and equating the 
coefficients at the like powers of the products e fe /i m 
in different combinations (A_ . . . A* _ . . .) to zero. In 
its essence, this procedure is similar to the method 
of undetermined coefficients. The coefficient calculation 
order and the general form of the expansions for A t 
and u n in arbitrary order of e are given in Appendix. 
The use of two parameters in ansatz lfl4|) is of key 
importance for the coefficient calculation procedure, 
since the expansions could not be split into linear- 
independent terms at e = /i. 

Note that our two-parameter approach has the same 
limitations in regard to convergence issues as other 
perturbative methods do. The convergence can get 
broken in the presence of resonances between harmonics, 
when quasi-harmonicity condition lfl3|) is violated. Some 
questions related to the convergence of asymptotic 
expansions for the solutions of differential equations were 
considered, in particular, in our works [18,19]. Expansion 
lfl4|) cannot be used either for wide-band wave trains 
with Afc ~ k. 

It should also be noted that the reduction of Eq. |(6|) 
with the second time derivative to Eq. I|14p with the first 
time derivative puts a constraint on the initial condition 
for u t . In this case, ut(x, 0) is a function of u(x, 0) defined 
by formula lfl4|) . 



4. High-Order Nonlinear Schrodinger Equation 

The two-parameter expansions were programmed in 
symbolic form for an arbitrary order of [i and e. The 



evolution equation for the complex amplitude of the first 
harmonic is 

A t = i.({±n)axA t + (ifi) 2 a 2 A^+ 

+(i^) 3 a 3 A^ + (i/x) 4 a 4 A^ + 0(a* 5 )+ 

+e 2 a 0j0i0 ^|A| 2 + (i/i) (ai.o.oA^Al 2 +a ,o,i^ 2 ^) + 

+ (in) 2 (a 2 ,o,oA^\A\ 2 +a hh0 A 2 A* + 

«o : o : 2A 2 A^ + 0(/i 3 ) + 
+e 4 a ,o,o,o,o^| 4 + 

+ (i/i) (ai,o,o,o,oAH^| 4 + a ,o,o,i,o^V| 2 ^) + 
+ (iM) 2 (a 2 ,o,o,o,o^«|^| 4 + ai,i,o,o,o^||^| 2 ^* + 
+a h o.o.i.olA^ 2 A\A\ 2 + a Q ^^ t2 . Q A 2 \A\ 2 A* a + 

+00,0,0,1,1^4) +°(^ 3 )] +°( £6 ))- (18) 

In each term of this equation, the power of the formal 
parameter [i points to the overall order of the derivatives 
with respect to £, and the power of the formal parameter 
e points to the nonlinearity order. These parameters 
disappear after going back to the original variables u\ = 
eA and x = Taking into account dispersion relation 
(fl2| . the coefficients a_ can be written as (at P = 5) 



a = 0, ai = , 



ft 2 



C 2 «i 



<13 = 



c A kot\ 
2lo 5 



0,4 



ai c 4 (4:c 2 k 2 - ax) 



8w 7 



ao,o,o 



Ol,0,0 



3«3 



5a 2 



2oj iuai ' 

_ c 2 k 
2fto. 0, 1 — — T 



1 d"^ 
n\ dfc" ' 



3a 3 



«2, 0, 



3ai 

— z(2c 2 k 2 (27a 1 a 3 - 14a: 



18a 

1 (27a 1 a 3 -62a 2 )), 



a i, 1,0 



36a 

ai(27aia3 — 38a^ 
2 



c 2 / 

=• 4c 2 fc 2 (27aia 3 - 28a 



a i,o, 1 



— - — = [Ac 2 k 2 (9aia 3 - 4a: 



-ai(9aia3 — 34a 2 , 
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c 2 / \ 
a °'°' 2 = 6a uj 5 ( c2fc2 ( 9aia3 + 2a 2) + lZctial); 

_ 1 / 335a| 25a| 143a^c 
a ,o,o,o,o - ~ 12, ,2^2 + 10 „,2 



+ 



w V 108a? 18u 2 af 12a? 
ha\a^ 9a 2 3a 2 , 14a 2 a4 



+ 



+ 



+ 5a 5 , 



2w 2 ai 8w 2 16ai ai 

c 2 fc/925a| 275a| 421ala 3 



01, 0,0, 0,0- ^3 Vl0ga 3 . ^2^2 ^2 

55a 2 a 3 99a 2 9a 2 42a 2 a 4 

2w 2 ai 8w 2 16ai ai 

c 2 fc/295a| lOOa^ 139a|a 3 



15a 5 ) : 



_ c"k f/ybaz lUUa 2 
20a|a 3 9a| 3a 2 , 28a 2 a 4 



u; 2 q;i 



6a 2 
10a 5 ) . 



8ai ai 

The expressions for the subsequent coefficients a_ are 
too long to be presented in explicit form. 

The complex amplitudes of other harmonics are 
found from the relations 

"o = £ 2 ([Cl^| 2 + + + 

+ ( iA1 ) 2 (6 2 0) A«^* + b^AsAl + b^ 2 AA^) + 



+0( M 3 )J + £ 2 [&o°o,o,ol^| 4 + 

+(i M ) (6^ 0i0 ^|A| 2 A* +^, A|A| 2 4) + 

+ (i/^) 2 (4°o, o, o% 1^* + &i°i, o, o^^* 2 + 

+ Cl,0lA| 2 |^| 2 + C2,0^l^| 2 ^ e + 

+Co,i,i^f)+°(^ 3 )] +0(e 4 )); 
"2 = e 2 ( [4V 2 + (iM)&i 2 U^+ 
+ (i Ai ) 2 (6 2 2 U«^ + &i 2 U|)+0(^ 3 ) _ 



+ 



+e 



b ( oio,oA 2 \A\ 2 + 
+(i M ) Klo.oA^I^I 2 + 4 2) o,o,i^|) + 
+ (i^) 2 (41 o, o^^l^l 2 + ^i, o, Al\A\ 2 + 
+b[ 2 ioA\ A ^ A2 + b oio,2A 3 A^) +0(» 3 )] +0(e 

«3 = £ 3 ([4 3 U^ 3 + (i^ 3) ,o^ 2 + 

+ (^) 2 (4 3) o,o^ 2 + &i!!o4M) + O^ 3 )} + 0(e 2 )) 
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At P = 5, the coefficients of the expansions 

(0) __202. (0) _ (0) 

°0.0 — ' °1.0 — °0, 1 — u : 

ai 



are 



ai 

h (0) _ 1 1,(0) _ do) _ 2c 2 a 2 
(o) 38a 2 10a 2 a 3 6a 4 

W,o = — ^3- + -^— -— ; 

JO) _ JO) _ . 
°1, 0,0,0 ~ °0,0, 1,0 — u > 

1,(0) _ & (0) 

u 2, 0, n n — u 

c 



'1, 0, 0, — °0, 0, 2, — 
.2 



c 2 / 

= 27q;3 ^ 4 (^-4c 2 £; 2 (23a 2 + 90aia 2 a 3 - 81a 2 a 4 ) + 
+a 1 (538a 3 , - 927a!a 2 a 3 + 324a?a 4 ) 



b (0) _ b (0) 

°1, 1,0,0 — "0,0, 1,1 — 

c 2 / 

= 27q,3^4 (~ 4c2fc2 ( 79a 2 + 18aia 2 a 3 - 81a 2 a 4 ) + 

+ai (314a^ - 639aia 2 a 3 + 324a 2 a 4 )), 
& (0) 

°1,0, 1,0 " 

4c 2 / 

= - a -. I -4c 2 fc 2 (l7a 2 + 18aia 2 a 3 - 27a?a 4 ) + 
9ajw 4 V 

+ai (I42a^ - 261aia 2 a 3 + 108a 2 a 4 )); 

,(2) _ _02_. ,(2) _ „. ,(2) _ _,(2) _ 2c 2 a 2 _ 
6 °> ~3ai' ^.°- U ' ° 2 ^- "Li" g^W 



( 2 ) _ 59a^ 31a 2 a 3 4a 4 

"0, 0, 0, — ~ ~ - - " I 



54a? 12a 2 3ai ' 

J 2 ) u( 2 ) 2c 2 fca 2 , 2 s 

6 i, o, o, o = & o, o, o, i = 27cj2q;3 (9ai"3 - 10a 2 ) ; 

2 

4 2) o,o,o = -j^s-ji (c 2 fc 2 (2606a 3 - 5283aia 2 a 3 + 
+1728a^a 4 ) +ai (5006af - 7443aia 2 a 3 + 1728a 2 a 4 )^) , 

2 

6 i 2)l '°'° = 129^7 (c 2 fc 2 (3574a 3 - 7191a 1 a 2 a 3 + 
+3456a^a 4 ) +ai (3094a^ - 6759aia 2 a 3 + 3456a 2 a 4 )^) , 

2 

Co,o,i = 162^4 (c 2 fc 2 (!22a 3 - 369a!a 2 a 3 + 
+432a 2 a 4 ) + a 4 (-358a 2 + 63aia 2 a 3 + 432a 2 a 4 )^) , 
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6 ( 2 ) - _E 

°0,0,0,2 - qo/1^,3,.,4 



(c 2 fc 2 ( 



2«2 — 261aiO!2C*3+ 



-432a 2 a 4 ) + ai (-238a^ - 45aia 2 a 3 + 432a 2 cV 



,(3) __^_,^3_. 6 (3) _ Q . 
3 o,o,o- 12a 2 + 8ai . fi.o.o-u, 



6 (3) _ _ b (3) 

°2, o, o — °i, 1, o 



(27aia 3 + 34a 2 .) 



288a 2 w 2 

Some of these coefficients were derived earlier in [9] 
by the multiple scale method. The expressions presented 
in [9] are in agreement with those obtained here (except 
for several misprints and typographic errors). 

5. The Effect of High-Order Dispersive Terms 

Let us illustrate the evolution of a wave train envelope 
described by the equation of form (f!4|) . To this end, we 
rewrite original equation (J6j) in dimensionless variables 
x = kx and t = ckt: 



p 



0; " = tv 



{ckf 



(19) 



In this case, we have c = 1, k = 1, up = ai + 1, and Uo 
is a typical amplitude of the function u. 

As an example, let us consider the case P = 3 with 
a.\ = 1, ct2 = 0, and S3 = —1/6. The values of these 
parameters correspond to the first two terms in the 
Taylor expansion of the function sin it. Hereafter, the 
tildes over the dimensionless variables are omitted. 

The corresponding coefficients of evolution equation 
QB are 



1 



1 



1 



"3 



8V2' 



ei4 



a o, o, o 



64^2' 
1 



) a o,o,i — 



8V5' 



0-2. o. o — ao. o. 2 



16\/2 



i a i, i,o — ai, o, l 



ao, o : o,o,o 



16^2' 

— V> ■■■ ■ ( 2 °) 
96^ 

Initially, we retain only those terms in Eq. QB whose 
overall order of smallness with respect to the parameters 
e and (i is no more than two. In this case, we obtain a 
classical NLSE: 



2|m| 

0.6- 

0.5 
0.4 
0.3 - 
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f = 



-4 
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4^ 
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" r r ' " T " 

100 200 300 400 500 600 700 800 900 1000 X 



Evolution of the wave train envelope which is given, at the initial 
moment t = 0, by function l|22|l with parameters (3 = ( = 1/10 and 
4>o = xo = 0. (Solid curve) exact one-soliton solution of NLSE (JHJ, 
(dashed curve) numerical solution of the generalized NLSE l|24|l 

It has an exact one-soliton solution at a 2 ao. 0.0 < 0: 

. 1/2 

u\(x, t) = (3 



\ a o, 0,0 



xexp ^i(-^=(x-ai<)-s(C 2 "/3 2 )t + 0o)J X 
x cosh -1 I Z— (x -x - ait) - 2s/?Ct ) . (22) 



02 



Here s = sign(ao,o,o) and (3, <fio, and xq are free 
parameters [20, 21]. The corresponding approximate 
solution of Eq. fl~9|) is 



(23) 



u(x, t) — ui(x, t) exp(i(\/2i — + 

+u{(x, t) exp(— i(\/2t — x)). 

To analyze the effect of high-order dispersive terms 
on the shape of one-soliton solution lj22]l , we consider the 
generalized NLSE 

(ui)t = -ai(ui) x - ±a 2 {ui) xx + a 3 (u 1 ) xxx + 
+ia 4 (Mi) 

xxxx i&O, 0, ^1 1^1 1 

-ai,o,o(wi)a:|wi| 2 - ao,o,iui{u*)x- 

-ia2,0,o(ui)a;£i;|ui| 2 - idi, i, ("l) 2 u\ — 

— iai,o, i|(wi)k| 2 wi - iao,o,2 u\(u\) xx + 



{u\) t = -ai{ui) x - ia 2 (ui) xx + ia ,o,oWiK| 2 . (21) + 



100,0,0,0,0 u l \ u l\ 



(24) 
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with the coefficients defined by (J20J) . The initial 
condition is chosen in the form of function l!22tl with 
/3 = £ = 1/10 and 4>o — xq — 0. The figure shows 
the evolution of such an envelope. To solve Eq. lf24"l) 
numerically, we used the split-step Fourier method [22, 
23]. High-order dispersive terms are seen to affect the 
amplitude, shape, and velocity of the soliton solution. 

6. Conclusion 

We described a general method for deriving the evolution 
equations for narrow-band wave trains in nonlinear 
media with dispersion. The procedure produces no 
secular terms and can easily be put in algorithmic form. 
By using the Klein-Gordon equation with arbitrary 
polynomial nonlinearity as an example, we derived a 
generalized NLSE whose coefficients can be calculated 
in any order with respect to nonlinearity and dispersion. 
The equation can be used to investigate the propagation 
of ultrashort pulses in optical fibers, to study wave events 
on a fluid surface, and to describe the Langmuir waves 
in hot plasmas. 

I.S. Gandzha thanks the State Fundamental 
Research Foundation of Ukraine for the financial support 
(grant of the President of Ukraine No. GP/F26/0056). 



papers devoted to high-order evolution equations. To our 
deep sorrow, these plans were ruined by the fatal malady 
and untimely death of Vasyl Petrovich. Only time will 
show whether the work which had started can be finished 
without its inspirer. 

Vasyl Petrovich was a kind, sincere, calm, even- 
tempered, tolerant, open, and generous man, who was 
always ready to help. He was faithful to his life principles 
and convictions till the end. He had ingenious non- 
standard way of thinking and well-trained intuition. 
He was a man of word and justice. Blessed memory 
about Vasyl Petrovich Lukomsky will always abide in 
the hearts of those people who had the honor to be in 
fellowship and collaboration with him. 

APPENDIX 

The general forms of the expansions for At and u n are 

oc oo 

n =0 k=0 

( M l)n , 2fc + l = 

no n\ n 2k — l 

= '■• "no-ni, "1-nj, ™2k-l-™2fc> n 2k X 

ni— 7i2~ n 2k = ® 
k+1 2k+l 

x I] A («i-i-m)i II ^K-j-^K- n 2k+1 =0; 



Postscript 

The major part of this paper was written during the 
last months of Vasyl Petrovich Lukomsky's lifetime 
(January 14, 1942 - March 31, 2008). The above- 
described two-parameter method was developed by him 
as far back as at the end of the 1980s (some of his 
notes regarding NLSE were dated by 1986 and 1987). 
Regretfully, the results of those studies were published 
with delay and only in one paper of 1995 [16]. Vasyl 
Petrovich was attracted by various fields of physics, he 
had very broad physical horizons and often sacrificed 
the time for paper preparation in favor of something 
new. Still unpublished is the paper "Method of indefinite 
coefficients for derivation of evolution equations with 
higher terms for nonlinear waves" written together with 
Yu.G. Rapoport and submitted to Physica Scripta in 
1999. They derived a high-order NLSE for describing 
the propagation of ion sound in non-isothermal plasmas. 
Vasyl Petrovich turned back to his method (which can 
appropriately be called the Lukomsky method) in 2005 
with the aim to derive a generalized NLSE for the 
modulations of gravity waves on a fluid surface. This 
paper was scheduled as the first out of the whole series of 



U n = £ n ^ (i/- 1 )" £2k ( U ")n ,2k + n , U > 2, 

n =0 k=0 
n ill 

(Un)n ,2k + n = X/ ' ' ' 

ni=0 n.2=0 

™2k + n-2 



■no, .... n 



2k + n-2 — n 2fc + Ti-l > n 2k + n-l 



X IT II A (n»-l-ni)«' n 2fc + n = 0; 

oo oo 

"0 = e 2 (i^)" e2fc ( u o)n , 2k+2 , 

n =0 k=0 

( u o)n ,2k + 2 = 

no "1 "2k 

= V V" V b {0) X 

l—i /—i "' Z—i "1— »2> —i "2fc— ™2fc + li n 2fc+l 

n 1= 0n 2 =0 n 2f! + i=0 

fc+1 2fc+2 

X IT II A (n«-l-™i)«> ™2fc + 2=0. 

i=l i = k + 2 

The sums above contain many identical summands (e.g., 
0,1,0, oA^AA* and oo, l, oAA% A*). The repetitions can be 
eliminated if one limits the summation orders using the rules 
n>i— 2 — ni-i ^ — rn ^ rij — ni+i. The corresponding 

expansions are 

(«l)n ,2fc + l = Yl Yl 

U\ =0 ri2= max(0, 2ni — no) 



ISSN 2071-0186. Ukr. J. Phys. 2009. V. 54, N 1-2 



213 



V.P. LUKOMSKY, I.S. GANDZHA 



Coefficient calculation order 



Iteration number 


Number of Eqs. (TTJ 


Coefficient 


k = (e 1 ) 


n = 1 (e 1 ) 










k = 1 (s 3 ) 


n = (e 2 ) 


6 (0) 

no— ni 


m 








n = 2 (e 2 ) 


ni 








n = 1 (e 3 ) 


^TlQ —Til 


ni— n 2 


n 2 




k = 2 (e 5 ) 


n = (e 4 ) 


fe (0) 


rti — n 2 


n 2 -n 3 , n 3 






n = 3 (e 3 ) 


b (3°r ni 

no— ni 


ni — n 2 


n 2 






n = 2 (e 4 ) 


ni — n 2 


n 2 -n 3 , n 3 






n = 1 (e 5 ) 




n 1 —n 2 


n 2 -n 3 , n 3 -n 4 


n 4 


k = 3 (e 7 ) 


n = (e 6 ) 


6 (o) 

n — ni 


ni —n 2 


n 2 -n 3 ,n 3 -ri 4 


ra 4 — 715 , 715 




n = 4 (e 4 ) 


ni —n 2 


n 2 -n 3 , n 3 






n = 3 (e 5 ) 


n\ —n 2 


n 2 -n 3 ,n 3 — n 4 


n 4 




n = 2 (e 6 ) 


n\ —n 2 


n 2 -n 3 ,n 3 — n 4 


n 4 — 715 ' n 5 




n = 1 (e 7 ) 


^TlQ — Tlx 


ni—n 2 


n 2 — n 3 , ?i 3 -n 4 


n 4 — n 5 , n 5 -n 6 , n 6 



E - E 

n-i =max(0, 2rij_ i — 2 ) n fc-f- 1 — max(0, 2n^. — i ) 



r 2fc + l-i ] 
[ 2fc + 2 -^ ^-l] 



E ••• E 

Tl fe + 2— Tij— max(0, 2n i _i — n i _ 2 ) 



[|«2fe-l] 

E 



'»0-»h ni-n 2 , ™2k-l-™2fc. ™2fc • 



••2* 



= max(0, 2n 2t _ 1 -n 2l _ 2 ) 



fc+1 2fc + l 

X EI E[ A K-l-7ii)€ ' «2fc + l=0; 



i=l 



i=fc + 2 



("n>l)n , 2fc + nlfc>0 — E E 

711=0 n 2 -max(0, 2ni — tiq) 

™i-l ™fc + n-l 

E •■■ E 

n i -niaxfO, 2nj_i — n j_ 2 ) n fc+ji- maJ( (0, 2n^_j_ Tl _i— 7ik + TZ _2 ) 



r 



[ afc '!tn-T+l "i-l] 



E ■•■ E 

ik + n + l ~ n i — max(0, 2rij_ ^ — _ 2 ) 

[^n 2fc + 71 _ 2 ] 2fc + n 

e *w n A _ n a : 

n 2fe + n-l= m ax(0, 2n 2fe + „_ 2 -"2fe + n-3) i = 1 i=fc + n+l 



L^-^ n 0j [ n _ i+ i »«-lJ 

(«n>l)n , 2fc + nU=0 = E "' E 

711=0 n j-max(0, 2n;_ i -nj_2 ) 

[i«-«-2] n 

e tW n^j 

n rx — 1 -maxfO, 2/1^ _ 2 — n 71 _3 ) i=l 
7i ni 

(«o)n , 2fc + 2 lfe>0 = E E 

711=0 ii 2 -max(0, 2ni — no) 

E •■• E 

rii =max(0, 2rii_i — n-i — 2) n k-\-l =max(0, 271/, — _ 1 ) 
r fc i r 2fc-i + 2 1 

E - E 

n fc + 2 — n-i — max(0, 2n^_ ^ — n^_ 2 ) 

[| n 2fc] + 1 2fc + 2 

e ^ ri A n 

™2fc + l= max (0, 27i 2fc -7i2fc_i) >=1 i=fe+2 

710 

("0)n ,2fc + 2lfc=0 = E ''no-Tii.ni^f'io-'ilJS^nie- 



The table gives the order of calculation of the coefficients &W 
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flBOnAPAMETPHHHHH METOfl flJIil OIHICy 
HEJIIHIHHOI EBOJIIOIIII CIIEKTPAJIbHO 
BY3BKHX XBHJIBOBHX I7AKETIB 



B.II. JlyKOMCbKUU 



I.C. Tandwca 



P e 3 io m e 



Po3rjiHHyTO iacoBy eBOjiioniio cneKTpajibHO By3bKiix xbhjibobhx 
naxeTiB KmneBOi aMnjiiTy^ii b HejiimiiHiii ^iicnepciftiiiii ciiCTe- 
Mi, onncyBamii piBHHHHHM KjieiSHa-LopflOHa 3 flOBijibHOio noni- 
HOMiajiLHOio HejiimiiHicTio. 3acTocoBaHO hobhh MeTO# TeopiT 36y- 
peHb, aKnii fl03BOjiae 3BecTH BnxiflHe xBiuibOBe piBHaHHa p,o mo- 
^ejibHoro piBHaHHa #jia o6Bi#HOi xBHjibOBoro naxeTy (HejiimiiHe 
piBHaHHa LUpefliHf epa BHinpro nopa^Ky). IIo6y r D;OBaHO acHMnTO- 
raiHe po3BHHeHHa lacoBOi noxiflHOi no abox He3ajie>KHHx napa- 
MeTpax, KOTpi xapaKTepH3yK>Tb MajiicTb aninjiiTyfl (e) i noBijib- 
HicTb ix 3MiH y npocTopi (fi). Ha BiflMmy Bifl iraniix MeTOfliB Teopi'i 
36ypeHb (TaKHx hk mcto^; 6araTbOx MacniTa6iB), #e ni ^Ba napa- 
MeTpn He po3pi3HaiOTb, flBonapaMeTpniHrrii MeTOfl He npnBOflHTb 
^o noaBH BiKOBHx (ceKyjiapHnx) ^o^aHKiB. Pe3yjibTaTH po6oth 
MO>KyTb 6yTH 3acTOCOBaHi #o ^ocjiifl^KeHHa nornnpeHHa yjiLTpa- 
KopoTKHx ((beMTOceKyHflHHx) iMnyjibciB b Mepeacax bojiokhhcto- 
onTHHHoro 3B'a3Ky, bhbighhh xBHjibOBHx aBHirj; Ha noBepxm pi^H- 
hh, onncy jieHrMiopiBCbKHX xbhjib b rapaiiii njia3Mi. 
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